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1.0  INTRODUCTION 


The  purpose  of  this  report  is  to  document  an  evaporation  duct  model  for  inclu¬ 
sion  as  a  Navy  standard  model  in  the  Oceanographic  and  Atmospheric  Master 
Library  (OAML). 

J  The  propagation  of  microwaves  at  altitudes  close  to  the  surface  of  the  sea 
depends  upon  the  refractive  structure  in  the  first  few  tens  of  meters.  This  height 
interval  is  called  the  atmospheric  surface  layer  and  is  usually  dominated  by  mechani¬ 
cal  and  buoyant  turbulence.  The  average  profile  of  wind,  temperature,  humidity,  or 
any  scalar  is  determined  by  the  turbulent  motions.  Roll'  presents  an  overview  of  the 
wind,  temperature,  and  moisture  fields  in  the  first  few  meters  above  t  he  sea  surface. 
Meteorological  measurements  to  directly  determine  these  profiles  in  the  surface  layer 
over  the  ocean  are  quite  difficult  and  essentially  impractical  on  an  operational  basis. 
To  circumvent  this  difficulty,  semiempirical  relationships  between  the  profiles  and 
the  fluxes  of  momentum,  heat,  and  moisture  have  been  developed.  , 

For  radio  purposes,  the  profile  of  potential  refractivity  in  the  surface  layer  is 
of  intercut,  with  potential  refractivity  defined  by 


77. h  /■„ 
0 


+  3.73 


eP 
0  1 


where  Pa  -  1000  mb,  o  is  potential  temperature  in  Kelvin,  and  ep  is  potential  vapor 
pressure  in  mb.  Further 


0  =  T 


(*r 


where  T  is  temperature  in  Kelvin,  P  is  pressure  in  mb,  and  e  is  vapor  pressure  in  mb. 
Over  the  ocean,  the  profile  of  potential  refractivity  is  greatly  shaped  by  the  presence 
of  trapping  gradients,  which  form  a  shallow  duct,  called  the  evaporation  duct  because 
of  the  predominant  influence  of  evaporation  in  causing  the  trapping  gradients.  The 
parameter  that  best  characterizes  propagation  is  the  evaporation  duct  height,  defined 
as  the  height  at  which  the  gradient  of  potential  refractivity  reaches  the  critical  value, 
calculated  by  ray  theory,  which  will  cause  a  ray  launched  horizontally  to  have  a  cur¬ 
vature  exactly  equal  to  the  curvature  of  the  earth.  In  terms  of  modified  refractivity, 
M,  the  evaporation  duct  height,  8,  is  the  height  in  the  surface  layer  at  which  the  M 
profile  reaches  its  minimum  value  and  dMUlz  =  0.  Figure  1  is  an  example  of  an  M 
profile  for  a  13-meter  evaporation  duct,  which  is  near  the  world  average  duct  height. 
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There  are  several  semiempirical  techniques  to  relate  surface  layer  profiles  to 
meteorological  measurements  that  can  be  simply  made  at  sea  (often  called  hulk  meas¬ 
urements).  Jeske  examined  several  schemes  and  correlated  evaporation  duct  height  to 
propagation  measurements.2"4  His  work,  with  modifications,5"'  is  the  basis  for  the 
duct  height  formulation  later  in  this  report;  this  formulation  is  developed  in  more 
detail  in  Appendixes  A,  B,  and  C.  A  FORTRAN  77  program  is  provided  in  Appendix 
D.  Rotheram8  has  developed  a  similar  formulation  in  a  slightly  different  manner. 


Figure  1.  Example  of  a  modified  reffactivity  profile  for 
a  13-meter  evaporation  duct.  At  the  surface  M=  325; 
z  is  the  altitude. 


2.0  INPUTS,  OUTPUTS,  AND  LIMITS 

The  bulk  meteorological  measurements  required  for  the  evaporation  duct 
calculation  are  wind  speed,  air  and  sea  temperature,  relative  humidity,  and  altitude 
of  the  measurements.  Since  gradients  in  the  surface  layer  generally  become  small  at 
altitudes  more  than  a  few  meters,  the  influence  of  measurement  altitude  is  neglected 
as  long  as  measurements  are  made  at  least  G  meters  above  the  surface.  Additionally, 
the  influence  of  sea  waves  on  the  measurements  is  usually  restricted  to  altitudes  less 
than  three  times  the  wave  height  (Ref.  6,  p.  123,  and  Ref.  9),  and  the  altitude  of  6 
meters  or  greater  minimizes  most  wave  effects. 

The  inputs,  outputs,  and  limits  are  discussed  in  the  following  paragraphs  and 
are  summarized  in  Table  1. 


Table  1.  Summary  of  inputs,  outputs,  and  limits. 


Parameter 

Type 

Limit 

Precision 

Wind  speed 

input 

0  to  50  kt 

1  kt 

Air  temperature 

input 

-20°  to  50°C 

0.1  °C 

Relative  humidity 

input 

0  to  100% 

1% 

Sea  temperature 

input 

0°  to  40°C 

0.1°C 

Evaporation  duct  height 

output 

0  to  40  ra 

0.1  ra 

2.1  Wind  Speed 

Wind  speed  can  be  measured  with  the  typical  cup  or  aerovaue  anemometer. 

The  anemometer  should  be  mounted  so  as  to  minimize  the  disturbance  to  air  flow 
caused  by  the  observing  platform.  On  a  ship  or  a  boat,  the  anemometer  should  be 
mounted  as  far  forward  and  as  high  as  possible;  the  measured  relative  wind  must  be 
corrected  for  ship’s  course  and  speed  to  determine  true  wind.  Kand-held  anemometers 
can  also  be  used  successfully  if  carried  to  the  windward  side  of  the  platform. 

The  wind  speed  should  be  averaged  over  a  period  of  5  to  10  minutes.  This  will 
average  out  the  turbulent  fluctuations  in  wind  speed.  Wind  speed  input  is  limited  to  0 
to  50  knots. 

2.2  Air  Temperature  and  Relative  Humidity 

Measurement  of  air  temperature  and  relative  humidity  can  be  made  at  any 
well-exposed  site  above  an  altitude  of  6  meters  with  reasonable  precautions  taken  to 
minimize  platform-induced  conductive  and  radiative  heating.  An  aspirated  psychro- 
meter  with  wet  and  dry  bulb  thermometers  is  a  convenient  portable  instrument. 

Since  temperature  and  humidity  gradients  above  6  meters  in  the  typical  overwater 
surface  layer  are  quite  small,  averaging  times  for  these  parameters  can  be  3  to  5 
minutes.  Air  temperature  input  is  limited  to  -20°  to  50°C;  relative  humidity  input  is 
limited  to  0  to  100%. 

2.3  Sea  Surface  Temperature 

The  simplest  measurement  of  sea  surface  temperature  is  made  with  a  small 
rubber  bucket  (or  any  container  with  a  small  thermal  mass)  and  an  accurate 
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thermometer.  More  durable,  specifically  designed  water  thermometer  buckets  are 
available  which  provide  a  frame  for  a  thermometer  and  a  small  reservoir  (bucket)  to 
retain  water  around  the  thermometer  bulb.  These  devices  are  used  to  sample  the 
water  over  a  finite  depth  (~  0.3  meter).  Under  most  conditions,  the  wind  is  sufficient 
to  mix  the  upper  levels  of  the  water  column  ( >  2  meters).  However,  on  cloudless  days 
with  strong  solar  heating  and  light  winds,  the  surface  temperature  can  vary 
considerably  from  the  underlying  water.  In  these  instances,  care  must  be  taken  to  get 
a  representative  sample  of  the  surface  water. 

Infrared  radiometer  sensors  provide  a  remote  means  to  sense  the  sea  surface 
temperature.  Although  considerably  more  expensive  than  the  buckets,  they  do  have 
the  advantage  of  sensing  only  the  sea  surface.  If  this  sort  of  instrument  is  used,  it 
should  have  a  means  to  compensate  for  reflected  sky  radiation.  Sea  temperature  input 
is  limited  to  0°  to  40°C. 

2.4  Evaporation  Duct  Height 

Based  on  the  inputs  of  wind  speed,  air  temperature,  relative  humidity,  and  sea 
surface  temperature,  the  evaporation  duct  height  is  calculated  as  the  only  output.  For 
purposes  of  propagation  assessment,  the  evaporation  duct  height  is  not  a  height 
below  which  an  antenna  must  be  located  to  give  enhanced  ranges  but  rather  a 
parameter  indicative  of  the  strength  of  the  evaporation  duct.  Strong  evaporation 
ducts  (high  duct  heights)  will  affect  surface-to-surface  propagation  at  frequencies  as 
low  as  1  to  2  GHz.  Weak  evaporation  ducts  flow  duct  heights)  will  affect  only  higher 
frequencies.  The  limits  on  evaporation  duct  height  are  0  to  40  meters. 
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3.0  MODEL 


The  following  six  steps  specify  the  calculation  of  evaporation  duct  height: 

1.  The  four  bulk  meteorological  measurements  required  ?re  <r  temperature, 
Ta  ,  and  seaj:  temperature,  Ts,  in  Celsius;  relative  humidity,  Rh,  u. .  -  •'■at;  and  wind 
speed,  u,  in  knots.  Temperatures  are  converted  to  Kelvin  by  Tat  -  T-  +-  273.2  and 
Tsk  =  Tj  +j  273.2.  If  wind  speed  is  less  than  0.01  knot,  evaporation  duct  height  is 

set  to  0  meters. 

;i 

2.  The  bulk  Richardson  number  is  calculated  by 
Rib  =j  369  2,  (Tak  -  Tsk)  /  (Tat  u1) 

where  z,  is  fhe  reference  height  in  meters.  Rib  is  restricted  to  be  no  greater  than  1. 

3.  Assuming  that,  in  the  surface  layer,  Tak  =  0  ,  e  -  ep,  and  neglecting  the 
effect  of  salinity  on  vapor  pressure,  calculate  potential  refractivity  for  the  air,  NA ,  at 
height  z,  and  for  the  air  immediately  in  contact  with  the  sea  surface,  N$: 


i 

77.6  j 

r 

1000  + 

4810 

Na  K 

Tak  1 

c 

Tak 

i 

r 

77.6  | 

r 

1000  + 

L 

4810 

Tsk  ' 

Tsk 

where  the  ambient  vapor  pressure  of  the  air,  e,  is  calculated  by 


!  RH 

e  =!  — —  es 
100 


es  =  6.105  exp  25.22 


Tak  -  273.2 


_  5.31  In  J-sO 
273.2 ) 


and  the  vapor  pressure  at  the  sea  surface,  e0 ,  is 
....  f _  Ttl  -  273.2 


e0  =  6.105  exp 


f  25.22  —  —  - 

I  T'* 


5.31  In 


4.  For  thermally  neutral  and  stable  conditions  (0  s  Rib  s  1),  the  evaporation 
duct  height.  6,  is  calculated  by 

(a)  8  =  0  for  A Np  a  0 

or 
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V>  = 

-4's  F 

F  *  -001 

v>  = 

10  ['  02  l0*  (^)  *  0  *’] 

-0.01  >  ■  2  -0.026 

v>  = 

10  [°-77s  l0*  (ir) 4  u  30‘3 

-0.026  >  J7  *  -0.1 

V  = 

10  10  630  lo*  (^0 4  0  “3 

1 

O 

*-+■ 

V 

^1- 

IV 

1 

V»  = 

io  t° 4,4  ,0*  (t1)  4  0  “3 

V 

IV 

1 

to 

to 

v>  = 

2 

7T  <  -2-2 
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6.  The  maximum  value  of  evaporation  duct  height  is  limited  to  40  meters. 
Remaining  variables  and  constants  are  defined  as  follows: 

ANp  =  Na  -  Ns  ,  potential  refractivity  difference  between  the  air  and  the 
sea 

-i  =  6  meters,  reference  height 

bi  =  -0.125  m  ' 1 ,  critical  gradient  of  potential  refractivity 
zo  =  1.5  x  10' 4  m  ,  surface  roughness  parameter 
g  =  4.5  ,  coefficient  in  the  unstable  function 

a  =  5.2  ,  coefficient  in  the  stable  function 

«  A-»  r 

l'  =  _ Lii  m ,  Monin-Obukhov  length,  where 

Rib 

r,  -=  0.05  Rib  s  -3.75 

r,  =  0.065  +  0.004  Rib  -3.75  <  Rib  £  -0.12 

T,  =  0.109  +  0.367  Rib  -0.12  <  Rib  S  0.14 

r,  ^  0.155  +  0.021  Rib  0.14  <  Rib 


The  dominant  factor  in  determining  evaporation  duct  height  is  the  difference 
in  potential  refractivity  between  the  air  at  the  reference  height  and  the  air  at  the  sea 
surface.  Errors  in  air  temperature  measurement  caused  by  conductive  and  radiative 
heating  effects  have  been  shown  to  strongly  affect  the  duct  height  calculation.7  To 
minimize  the  sensitivity  of  the  evaporation  duct  algorithm,  an  additional  test  is 
applied  whenever  T0  -  Ts  >  -  1  .  Evaporation  duct  height  is  calculated  for  Ta  =  Ts , 
(d0),  and  for  Ta  =  T,  -  1,  with  Ts ,  u,  and  RH  unchanged.  Then,  if  60  >  d.u  the 
value  of  <J_i  is  the  evaporation  duct  height;  otherwise,  the  evaporation  duct  height  is 
calculated  as  defined  in  steps  1  through  6  above. 
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4.0  TEST  CASES  FOR  CALCULATION  OF 
EVAPORATION  DUCT  HEIGHT 

Table  2  lists  data  that  can  be  used  to  verify  the  proper  implementation  of  the 
computer  program  in  Appendix  D.  For  the  specified  inputs  of  air  temperature  (°C), 
sea  temperature  (°C),  relative  humidity  (%),  and  wind  speed  (knots),  the  evaporation 
duct  should  be  calculated  correctly  to  within  ±.0.1  meter. 


Table  2.  Test  inputs  and  evaporation  duct  height  (8)  outputs. 


Ta 

(°C) 

— 

RH 

(%) 

U 

(knots) 

6 

(m) 

0.0 

15.0 

10 

1 

3.9 

0.0 

15.0 

20 

16.8 

0.0 

15.0 

75 

1 

2.8 

0.0 

15.0 

75 

20 

12.3 

0.0 

15.0 

100 

1 

2.4 

0.0 

15.0 

100 

20 

10.4 

0.0 

16.0 

10 

1 

4.0 

0.0 

16.0 

10 

20 

17.3 

0.0 

16.0 

75 

1 

3.0 

0.0 

16.0 

75 

20 

13.0 

0.0 

16.0 

100 

1 

2.6 

0.0 

16.0 

100 

20 

11.2 

15.0 

15.0 

10 

1 

8.7 

15.0 

15.0 

10 

20 

37.0 

15.0 

15.0 

75 

1 

3.6 

15.0 

15.0 

75 

20 

13.9 

15.0 

15.0 

100 

1 

0.0 

15.0 

15.0 

100 

20 

0.0 

15.0 

16.0 

10 

1 

9.1 

15.0 

16.0 

10 

20 

38.8 

15.0 

16.0 

75 

0 

0.0 

15.0 

16.0 

75 

1 

3.8 

15.0 

16.0 

75 

20 

14.6 

15.0 

16.0 

100 

1 

0.9 

15.0 

16.0 

100 

20 

2.6 

30.0 

15.0 

10 

1 

8.7 

30.0 

15.0 

10 

20 

37.0 

30.0 

15.0 

75 

1 

3.6 

30.0 

15.0 

75 

20 

13.9 

30.0 

15.0 

100 

1 

0.0 

30.0 

15.0 

100 

20 

0.0 

30.0 

16.0 

10 

1 

9.1 

30.0 

16.0 

10 

20 

38.8 

30.0 

16.0 

75 

1 

3.8 

30.0 

16.0 

75 

20 

14.6 

30.0 

16.0 

100 

1 

0.0 

30.0 

16.0 

100 

20 

0.0 
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Appendix  A 


DEVELOPMENT  OF  THE  EVAPORATION  DUCT 
HEIGHT  FORMULATION 

References  A-l  to  A-3  provide  a  detailed  development  of  the  evaporation  duct 
height  formulation  by  Jeske.  A  general  development  of  Jeske’s  formulation  is 
included  here  along  with  assumptions  and  rationale  for  the  implementation  in  the 
computer  program  listing  in  Appendix  D. 

Jeske  defined  the  gradient 


d* 

dz 


s, 

Q  <  U-  (z  +  ZQ) 


<t> 


(A-l) 


where  i  is  a  scalar  property  of  the  atmosphere,  S,  is  the  vertical  flux  of  t ,  z  is  height 
in  meters,  k  is  von  Karmen’s  constant  (0.4),  u •  is  the  friction  velocity  in  m/s,  za  is  the 
roughness  parameter  in  meters,  and  0  is  a  stability-dependent  function  (equal  to  1  at 
thermally  neutral  conditions,  i.e.,  air  temperature  equal  to  sea  surface  temperature). 
For  stable  conditions  (air  warmer  than  water),  the  stability  function  was  taken  to  be 
the  logarithmic-linear  model  proposed  by  Monin  and  Obukhov 


(A-2) 


where  a  is  taken  to  be  5.2  and  L'  is  the  gradient  form  of  the  Monin-Obukhov  scaling 
length  corrected  for  stability.  For  unstable  conditions  (water  warmer  than  air),  the 
KEYPS  relationship 


<P*  ~ 


=  1 


(A-3) 


with  a  -  4.0  is  used.  To  evaluate  z/L' ,  the  bulk  Richardson  number  is  used  as 


Rib 


X*Tz 
T  u2  r 


(A-4) 


where  g  =  9.8  m/s2,  z  is  height  in  meters,  T  is  air  temperature  in  Kelvin,  AT  is  air 
temperature  minus  sea  temperature,  u  is  wind  speed  in  m/s,  and  T  is  the  profile  coef¬ 
ficient.  Rib  involves  approximating  gradients  with  differences  and  assumes  T  -  6  in 
the  surface  layer.  The  profile  coefficient  for  neutral  conditions  is 


2i  ♦  Zq\  (A-5) 

2"  J 


A-l 


l 


and  is  approximately  0.1  for  typical  observation  heights,  zi,  and  surface  roughness 
parameter  values,  :a,  of  0.00015  meter.  Similarly,  L'  is 


V 


r  u 2  re 

g  AT 


(A-6) 


or,  in  terms  of  Rib 


_1_  _  Rib 
L'  '  io  z,  r. 


(A-7) 


where  re  is  an  empirical  profile  coefficient  for  nonneutral  conditions  (Fig.  A-l).  F,  is 
evaluated  by  a  straight-line-segment  fit  to  the  data  in  Fig.  A-l. 


3  00  2  SO  2  00  - 1  50  - 1  00  -0  50  0  00  0  50  1  00 

Rib 

Figure  A-1 .  The  profile  coefficient,  r, ,  as  a  function  of  bulk  Richardson 
number,  mb  (from  JeskelA_1  through  A-3)  j 


Returning  to  Eq.  A-l,  assume  the  flux,  S, ,  is  constant  with  height  and  inte¬ 
grate  from  the  sea  surface  to  z\. 


^  S‘  f  *  A 

ii  -  i0  =  An  =  - - dz 

Q  K  U-  J  Z  *Z0 
0 


(A-8) 


A-2 


Thus,  the  flux  term  is 


S, 

g  k  u* 


Aii 


I 


Z  +  2a 


and  substituting  in  Eq.  A-l  yields 


di  At  0 

dz  B(z  +  z0) 

where 


B 


/ 


Z+20 


(A-9) 


(A- 10) 


(A-ll) 


Taking  the  scalar  <  as  potential  refractivity,  Np ,  the  height,  6,  at  which  the  critical 
gradient  of  potential  refractivity  (bx  =  -0.125  m_1)  occurs  is  found  from  Eq.  A-10  as 


6 


ANp 

b\B  -A Np~ 


(A-12) 


where  <P  for  stable  conditions  is  determined  from  Eq.  A-2  and  B  is  evaluated  as 


B 


(A- 13) 


Jeske  recommends  that  6/L'  should  not  exceed  1.  If  this  occurs,  then  duct  height 
should  be  recalculated  with  6/L'  -  1. 


For  unstable  conditions,  solving  Eq.  A-10  for  0  in  terms  of  Np,  substituting  in 
Eq.  A-3  with  z  =  6  ,  and  solving  yields 


f  b'B\  i£  f mX " 7 

"  L'  ^A/VpJ 


(A- 14) 


where  B  is  evaluated  as 


V<Z|/L')  has  the  analytic  form*-4 


(A- 15) 


V>  =  1-^-3  ln0  +  2  ln| 


0?) 


+  2  tan'1  0  - 


n 

~2 


(A-1G) 


A-3 


Use  of  Eq.  A- 16  requires  the  iterative  solution  of  Eq.  A-3  for  ^  with  a  given 
z/L’.  'P  can  be  more  quickly  determined  by  a  straight-line-segment  fit  to  theKEYPS 
profile  of  Fig.  A-2  adapted  from  Lumley  and  Panofsky.A*5 


Figure  A-2.  Universal  function  *  as  a  function  at  z/L'  (from 
Lumley  and  Panofsky  A ' 5 ). 
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Appendix  B 


CRITICAL  GRADIENT  OF  POTENTIAL  REFRACTIVITY 

The  equation  for  refractivity  in  terms  of  pressure,  P,  in  mb,  temperature,  T,  in 
Kelvin,  and  water  vapor  pressure,  e,  in  mb  is 

N  =  ~ -  +  3.73  x  105  -L-  (B-1) 

For  application  in  the  atmospheric  surface  layer,  potential  refractivity 

Np  =  *  3.73  x  10s  ^  (B-2) 

is  the  conservative  property  where  P0  is  a  reference  pressure  level  (taken  to  be  1C00 
mb),  e  is  potential  temperature,  and  ep  is  potential  water  vapor  pressure. 

From  geometric  optics,  the  critical  gradient  required  for  trapping  is  that 
which  yields  a  ray  curvature  equal  to  the  earth’s  curvature 


—  - - =  -0.157  Nlm  (B-3) 

dz  a 

where  a  is  the  earth's  radius,  6371  x  103  meters.  To  relate  the  gradient  of  potential 
refractivity  to  the  gradient  of  refractivity,  take  the  total  derivatives  of  Np{Q.  ep)  and 
N(P,  T,  e)  with  respect  to  z.  Assuming  0  ~  t  and  ep~t  in  the  surface  layer  yields 

dNp  dN  dN  dP 

dz  dz  dP  dz 

The  partial  derivative  dN/dP  =  77.6 IT  varies  from  0.28  to  0.26  mb"1  over  the  range  of 
temperature  0°  to  30°C.  The  derivative  dP/dz  can  be  evaluated  by  making  the  hydro¬ 
static  approximation 


=  -eg  =  -- 


where  q  is  density,  g  is  acceleration  of  gravity  (9.8  m/s2),  P  is  pressure  (1000  mb),  R 
is  the  individual  gas  constant  for  dry  air  (2.87  x  106  erg/g/°K),  and  T  is  temperature 
(°K).  This  yields  a  variation  of  dP/dz  from  -0.12  to  -0.11  mb/m  over  the  temperature 
range  of  0°  to  30°C.  Taking  standard  temperature  (15°C/,  me  critical  gradient  for 
trapping,  in  terms  of  potential  refractivity  is 


-£■  =  -0.157  -  (0 . 27) (—  0. 12)  =  - 0. 125m'1 
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Appendix  C 


IMPLEMENTING  A  MODIFIED  EVAPORATION  DUCT 
HEIGHT  CALCULATION 

Reference  C-l  showed  that  a  bias  in  air-sea  temperature  difference  toward 
thermal  stability  resulted  in  spuriously  high  calculated  evaporation  duct  heights.  A 
modification  to  the  duct  height  calculation  that  detected  and  compensated  for  the 
air-sea  temperature  difference  bias  was  proposed  based  on  the  variation  of  duct 
height  as  a  function  of  air-sea  temperature  difference  at  a  constant  relative  humidity. 
To  implement  the  algorithm,  the  value  of  either  the  air  temperature  or  sea  tempera¬ 
ture  may  have  to  be  changed  to  yield  a  modified  air-sea  temperature  difference.  A 
comparison  of  climatological  air  and  sea  temperatures  from  Marsden  square  data0  - 2 
and  ocean  data  buoys0- 3  showed  better  agreement  between  sea  temperature  means 
than  between  air  temperature  means  (Tables  C-l  through  C-5).  Air  temperature 
means  from  the  Marsden  square  data  also  indicate  a  diurnal  variation  (~1°C)  that 
the  buoy  data  do  not  show.  Mean  air  temperatures  for  3-hour  intervals  for  data  buoys 
show  diurnal  variations  of  only  a  few  tenths  of  a  degree.*  This  gives  more  weight  to 
holding  sea  temperature  constant  and  varying  air  temperature  to  obtain  the  appropri¬ 
ate  air-sea  temperature  difference. 

The  next  question  that  arises,  then,  is  whether  or  not  it  is  reasonable  to  hold 
relative  humidity  constant  while  changing  air  temperature.  Reed  reported  finding 
systematic  errors  in  air  temperature  measurements  aboard  a  ship  but  did  not  find 
systematic  errors  in  relative  humidity  measurements.0- 4  Experience  with  psychro- 
metric  measurements  in  a  thermal  screen  during  the  Ku-Band  Surface  Surveillance 
Project0  5  sited  on  Point  Loma  in  San  Diego,  tends  to  support  Reed’s  findings.  Rec¬ 
ommendations  of  a  working  group  on  the  computation  of  global  air-sea  flux  climatol¬ 
ogy  to  use  both  daytime  and  nighttime  humidities  but  only  nighttime  air  and  sea 
temperature  lend  support  to  the  assumption  of  negligible  bias  in  the  relative  humid¬ 
ity  measurements.0  - 4 

The  implementation  of  the  modified  evaporation  duct  height  calculation 
assumes  that  air  temperature  can  be  varied  to  obtain  the  appropriate  air-sea  tem¬ 
perature  differences  while  sea  temperature,  wind  speed,  and  relative  humidity  are 
held  constant.  Figure  C-l  is  a  flow  chart  of  the  process. 


*  D.  B.  Gilhousen,  National  Oceanic  und  Atmospheric  Administration  Data  Buoy 
Center,  personal  communication,  January  1984. 
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Table  C-l.  Mean  air  and  sea  temperatures  in  °C  for  Marsden  square  81  (day 
and  night)  and  NOAA  data  buoy  42003  located  26.0  N  86.0  W. 


Air 

Sea 

Month 

MS81  ] 

42003 

MS81 

|  42003 

Jan 

22.2 

19.4 

24.0 

23.7 

21.2 

24.0 

Feb 

22.0 

19.8 

23.6 

24.2 

20.7 

23.5 

Mar 

23.8 

21.9 

24.1 

24.8 

22.4 

23.9 

Apr 

24.7 

23.9 

25.1 

25.7 

23.7 

25.0 

May 

26.5 

25.4 

26.6 

26.6 

25.5 

26.5 

Jun 

28.0 

27.7 

28.2 

?8.5 

27.0 

28.2 

Jul 

28.9 

28.2 

29.1 

29.3 

27.9 

29.1 

Aug 

28.9 

28.2 

29.4 

29.3 

28.0 

29.4 

Sep 

28.6 

28.0 

29.0 

29.1 

27.6 

29.0 

Oct 

27.1 

25.9 

27.9 

27.9 

26.0 

27.9 

Nov 

25.0 

23.7 

26.4 

27.0 

23.9 

26.4 

Dec 

23.4 

21.5 

25.1 

25.9 

22.3 

25.1 

C-2 


\ 


Air 


Sea 


Month 


MS194 


46004 


MS194 


46004 


Table  C-4  Mean  air  and  sea  temperatures  in  °C  for  Marsden  square  195  (day 
and  night)  and  NOAA  data  buoy  46001  located  56.0  N  148.0  W. 


Air 

Sea 

Month 

MS195 

46001  j 

MS195 

46001 

Jan 

4.6 

2.6 

5.4 

4.0 

4.2 

5.3 

Feb 

4.2 

2.2 

4.9 

3.6 

3.5 

4.8 

Mar 

4.7 

2.5 

4.8 

3.7 

3.6 

4.7 

Apr 

5.5 

3.6 

5.0 

4.6 

4.7 

5.2 

May 

7.2 

4.9 

6.1 

5.6 

6.1 

6.1 

Jun 

9.5 

8.3 

8.0 

9.1 

8.2 

7.8 

Jul 

12.0 

11.0 

10.6 

11.7 

10.7 

10.4 

Aug 

13.4 

12.1 

12.3 

12.7 

12.2 

12.4 

Sep 

12.7 

10.9 

12.2 

11.6 

11.5 

11.9 

Oct 

9.4 

7.6 

9.7 

8.6 

8.5 

9.4 

Nov 

6.5 

5.0 

7.6 

6.6 

5.9 

7.3 

Dec 

4.8 

2.4 

6.2 

4.9 

4.3 

6.0 

Table  C-5.  Mean  air  and  sea  temperatures  in  °C  for  Marsden  square  196  (day 
and  night)  and  NOAA  data  buoy  46003  located  52.0  N  156.0  W. 


Air 

Sea 

Month 

MS196 

46003 

MS196 

|  46003 

Jan 

3.8 

3.5 

4.4 

4.4 

3.4 

4.4 

Feh 

3.2 

1.8 

3.9 

3.6 

2.3 

3.7 

Mar 

3.6 

3.0 

3.8 

3.5 

2.7 

3.7 

Apr 

4.5 

3.7 

4.0 

4.1 

3.3 

3.9 

May 

6.4 

4.7 

5.0 

5.1 

4.9 

4.9 

Jun 

8.7 

7.2 

6.8 

7.3 

7.1 

6.6 

Jul 

11.3 

9.« 

9.6 

10.1 

9.8 

9.4 

Aug 

12.6 

11.2 

11.3 

11.5 

11.3 

11.1 

Sep 

11.9 

10.7 

11.1 

11.5 

10.5 

10.9 

Oct 

8.7 

7.6 

9.0 

9.0 

7.7 

8.9 

Nov 

5.8 

4.8 

6.6 

6.6 

5.0 

6.7 

Dec 

4.2 

3.2 

52 

5.0 

3.6 

0.1 
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Appendix  D 


FORTRAN  PROGRAM  TO  CALCULATE  EVAPORATION 

DUCT  HEIGHT 

c 

c  **********************<■******************************* 

c  name:  delta.f 

c  routine  controls  the  execution  of  the  ’get_delta’ 
c  function  for  testing  the  calculations  of  the  evaporation 
c  duct  height  (delta). 

c  language:  fortran  iv  (’77) 

c  machine:  hp  520  (unix) 

c  call  sequence: 

c  M4IN  routine 

c  glossaiy: 

c  airt:  air  temperature,  deg  C 

c  seat:  sea  temperature,  deg  C 

c  rh:  relative  humidity,  % 

c  u:  wind  speed,  knots 

c  subroutines: 
c  get  delta  (function) 

c  evap 

c  rev  date  description 

c  0.00  101987 

c  ***+********#************•«*******************«******» 

c  prompt  the  user  for  entries,  call  'get_delta'  to 
c  compute  the  evaporation  duct  height 

c  ******«**••••*•*******«*••*•**•*•****•*•*«****•***«•** 

write!*, 9000) 

9000  format(  “Compute  the  evaporation  duct  height  from"/ 

*  “entries  of:”/ 

*  “  air  temperature,  deg  C7 

*  “  sea  temperature,  deg  C7 

*  “  relative  humidity,  %  and"/ 

*  “  wind  speed,  knots’/) 


D-l 


/ 


write(7,9010) 

9010  format(“.  .  air. .  .|.  sea  . .  .|. .  rh. .  ,|. .  .  u  . . .  | .delta.) 7 
*  “  temp  |  temp  )  %  |  knots  |  meters") 


do  while(.true.) 

write{*, ’(“enter  air  temp  (C):  ",$)’) 
read(*,*)  airt 

write(*, ’(“enter  sea  temp  (C):  ”,$)’) 
read(*,*)  seat 

write{*, ’(“enter  relative  humidity  (%):  *,$)’) 
read(*,*)  rh 

write(*, ’(“enter  wind  speed  (knots):  *,$)’) 
read(*,*)  u 

delta =get_delta(airt,  seat,  rh,  u) 
write(7,9010) 

write(7,9020)  airt, seat, rh,u, delta 
9020  format(5f8.1) 

enddo 


end 


******«*********************************************** 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


name:  get_delta.f 

routine  calculates  the  evaporation  duct  height  using 
the  paulus  formulation  "Practical  Application  of  the 
IREPS  Evaporation  Duct  Model",  NOSC  TR  966,  Jun  84 

language:  fortran  iv  (’77) 
machine:  hp  520  (unix) 

call  sequence: 

delta  =  get_delta(airt,  seat,  rh,  u) 
giossaiy: 

airt:  air  temperature,  deg  C 

seat:  sea  temperature,  deg  C 

rh:  relative  humidity,  % 

u:  wind  speed,  knots 

subroutines: 

evap 

rev  date  description 

0.00  101987  modify  ’duct63’  program  code  to 
make  a  stand-alone  utility. 


function  get_delta(airt,seat,rh,u) 


c  check  wind  speed,  convert  to  deg  K,  and  compute 
c  vapor  pressures 

if  (u  .gt.  .01)  then 

wvel  =  amax  1(1.0,  u) 

dtok  =  273.2 

tak  =  airt  +  dtok 
tsk  =  seat  +  dtok 

esw  *  6. 105*exp(25.22*seat'tsk-5.31*alog(tsk/dtok)) 
es  ■  6. 105*exp(25.22*airt/tak-5.3l*alog(takAltok)) 
e  «rh*es/100. 


D-3 


c  calculate  delta 


tl  =  tsk-1.0  [modified  air  tem 

ifitak  .gt.  tl)  then 
el  =  rh*esw/100. 

call  evap(wvel,tsk,tsk,el,e8w,del0) 

esl  =  6.105*exp(25.22*(seat-1.0)/tl-5.31*alog(tl/dtok)) 

el  =  rh*esl/100. 

call  evap(wvel,tl,tak,el,esw, delta) 
ifidelO-delta  .It.  0.)  then 
call  evap(wvel,tak,tak,e,esw, delta) 
endif 
else 

call  evap(wvel,tak,tsk,e,esw,della) 
endif 


delta  =  0. 
end  if 
return 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


***************************************************** 

Name:  EVAP 

Routine  calculates  the  evaporation  duct  height 
(in  metres)  following  the  method  of  Jeske  "The  State 
of  Radar  Range  Prediction  over  Sea",  Tropospheric 
Radio  Wave  Propagation  -  Part  II,  NATO-AOARD,  Feb  71 

Language:  FORTRAN  IV  (’77) 

Machine :  Univac  1100/82 

CALL  Evap(wvel,tak,tsk,e, esw, delta) 

wvel:  Wind  speed  knots 

tak:  Air  temperature  in  Kelvin 

tsk  :  Sea  temperature  in  Kelvin 

e  :  Vapor  pressure  mbs 

esw  :  Surface  vapor  pressure  mbs 

delta  Duct  height  m  RETURNED 

Rev  Date  Description 
0.0  052984 


subroutine  evap(wvel, tak, tsk,e, esw  .delta) 

real lnzlzO 

delta  =  0. 

{default  duct  height 

zO  =  1.5e-4 

{roughness  height 

zl  =  6. 

{reference  height  (m) 

beta  =  4.5 

{constant 

alpha  =  5.2 

{constant 

lnzlzO  =  alog(zl/z0) 

{useful 

bl  =  -0.125 

{potential  refractivity  trapping  grad. 

rib = 369 .  *  z  1  *  ( tak-tsk)/(tak*  wvel  *  wvel) 
iftrib  .gt.  1.)  rib  =  l. 
c 

c  compute  potential  refractivity  difference 
c 

delna= 77.6/tak*(1000. + 4810.Vtak) 
delnO  =  77.6/tsk*(1000. + 4810.*esw/tsk) 
deln  =  delna-delnO 
iffdeln  .ge.  0.)  return 
c 

c  - 

c  compute  the  gamma  function 
c 

iffrib  .le.  -3.75)  then 
gamma =0.05 
else  iffrib  .le.  -0.12)  then 
gamma = 0.065  +  rib*0.004 
else  iffrib  .le.  0.14)  then 
gamma = 0.109  +  rib*0.367 
else 

gamma = 0. 155  +  rib*0.02 1 
end  if 
c 

olp  =  rib/(10.*zl*gamma)  !1./L' 

zlolp  =  zl*olp 
c 

c - 

c  compute  Psi  function  for  unstable  conditions 
c 

iffrib  .It.  0.)  then 

iffzlolp  .ge.  -0.010)  then 
psi=-4.5*zlolp 
else  iffzlolp  .ge.  -0.026)  then 

psi  =  10.**(1.020*alogl0(-zlolp)  +  0.690) 
else  iffzlolp  .ge.  -0.100)  then 

psi  -  10.**(0.776*alogl0(-zlolp)  +  0.306) 
else  iffzlolp  .ge.  -1.000)  then 

psi  =  10.**(0.630*alogl0(-zlolp) +0.160) 
else  iffzlolp  .ge.  -2.200)  then 

psi*  10.**(0.414*alogl0(-zlolp)  + 0.160) 
else 

psi  =  2.000 
end  if 
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c 

c  Compute  unstable  delta 
c 

b  =  lnzlzO-psi 
btemp  =  b*bl/deln 

d  =  btemp**4  -  4.*beta*olp*btemp**3 
ifXd  .gt.  0.)  delta =d**(-0.25) 
c 

c  Stable  and  neutral  conditions 
c 

else 

b  =  lnzlzO  +  zlolp*alpha 

delta  =  deln/(b*bl  -  alpha*deln*olp) 

iftdelta  .It.  0.  .or.  delta*olp  .gt.  1.)  then 

delta  =  (deln*(l .  +  alpha)-b  1  *alpha*z  l)/(b  1  *  Inz  1  zO) 
end  if 
end  if 
c 

c  Return  to  caller 
c 

return 

c 

end 
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